AFRL-PR-WP-TP-2007-222 


COOLING  PERFORMANCE  OF  A 
PARTIALLY-CONFINED  FC-72 
SPRAY:  THE  EFFECT  OF 
VARIABLE  GRAVITY  (Postprint) 


Travis  E.  Michalak,  Kirk  L.  Yerkes,  Rebekah  Puterbaugh, 
Scott  K.  Thomas,  and  John  McQuillen 


JANUARY  2007 


Approved  for  public  release;  distribution  is  unlimited. 


STINFO  COPY 


The  U.S.  Government  is  joint  author  of  the  work  and  has  the  right  to  use,  modify,  reproduce, 
release,  perform,  display,  or  disclose  the  work. 


PROPULSION  DIRECTORATE 

AIR  FORCE  RESEARCH  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7251 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 

Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YY) 

January  2007 

2.  REPORT  TYPE 

Conference  Paper  Postprint 

3.  DATES  COVERED  (From  -  To) 

07/01/2004-12/15/2006 

4.  TITLE  AND  SUBTITLE 

COOLING  PERFORMANCE  OF  A  PARTIALLY-CONFINED  FC-72 

SPRAY:  THE  EFFECT  OF  VARIABLE  GRAVITY  (Postprint) 

5a.  CONTRACT  NUMBER 

In-House 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

62203F 

6.  AUTHOR(S) 

Travis  E.  Michalak,  Kirk  L.  Yerkes,  and  Rebekah  Puterbaugh  (Power  Division 
(AFRL/PRP)) 

Scott  K.  Thomas  (Wright  State  University) 

John  McQuillen  (NASA  Glenn  Research  Center) 

5d.  PROJECT  NUMBER 

3145 

5e.  TASK  NUMBER 

20 

5f.  WORK  UNIT  NUMBER 

314520C9 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Power  Division  (AFRL/PRP) 

Propulsion  Directorate 

Air  Force  Research  Laboratory,  Air  Force  Materiel  Command 
Wright-Patterson  AFB,  OH  45433-7251 

Wright  State  University 

Dayton,  OH  45435 

NASA  Glenn  Research  Center 
Cleveland,  OH  44135 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

AFRL-PR-WP-TP-2007-222 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Propulsion  Directorate 

Air  Force  Research  Laboratory 

Air  Force  Materiel  Command 

Wright-Patterson  AFB,  OH  45433-7251 

10.  SPONSORING/MONITORING  AGENCY 
ACRONYM(S) 

AFRL-PR-WP 

11.  SPONSORING/MONITORING  AGENCY 
REPORT  NUMBER(S) 

AFRL-PR-WP-TP-2007-222 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

Presented  at  the  45th  AIAA  Aerospace  Sciences  Meeting  and  Exhibit,  Grand  Sierra  Resort  Hotel,  Reno,  NV,  8-11  January  2007. 
PAO  Case  Number:  AFRL/WS  06-2859,  14  December  2006.  The  U.S.  Government  is  joint  author  of  the  work  and  has  the  right  to 
use,  modify,  reproduce,  release,  perform,  display,  or  disclose  the  work. _ 

14.  ABSTRACT 

This  paper  discusses  the  effects  of  a  variable-gravity  environment  on  the  performance  of  a  subcooled  partially-confined 
spray.  An  experiment  was  flown  on  the  NASA  KC-135,  providing  various  gravity  levels  (-  0.16,  0.35,1.0,  and  1.8  g). 
The  experiment  consists  of  a  nozzle  spraying  on  a  Thick  Film  Resistor  (TFR)  heater  which  is  mounted  on  an  insulating 
glass  pedestal.  The  working  fluid  is  FC-n.  Non-dimensional  heat  input  varied  over  the  range  35  <  (1-f)  GDelta  <  190, 
with  a  heat  flux  of  20  <  (l-f)q"  <  70  W/cm2.  The  non-dimensional  grouping  (Frl/2Ga)l/2,  a  ratio  of  inertia  and 
acceleration  forces  to  viscous  forces,  varied  over  the  range  35  <  (Frl/2Ga)l/2  <  70.  Subcooling  temperature  ranged  from 
24°  C  <  Tsc  <  30°C.  For  this  data,  from  one  flight  week  of  three  flights,  DeltaT  =  Ts-Tsat  is  found  to  decrease  with 
decreasing  values  of  (Frl/2Ga)l/2.  It  appears  that  higher  values  for  subcooling  may  enhance  this  acceleration  effect. 
Changing  the  flow  rate  may  affect  how  much  the  acceleration  level  affects  the  2-phase  heat  transfer  coefficient  (h2-p), 
with  the  higher  flow  rate  showing  less  of  a  change  in  h2-p  as  the  acceleration  level  changes. _ 


15.  SUBJECT  TERMS 

Thermal  Management,  Spray  Cooling,  Electronics  Cooling,  Variable  Gravity 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT: 

SAR 

18.  NUMBER 

OF  PAGES 

20 

19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

Kirk  L.  Yerkes 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

(937)  255-6226 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


1 


45th  AIAA  Aerospace  Sciences  Meeting  and  Exhibit 
8-11  January  2007,  Reno,  Nevada 


AIAA  2007-198 


Cooling  Performance  of  a  Partially-Confined  FC-72  Spray: 

The  Effect  of  Variable  Gravity 

Travis  E.  Michalak*,  Kirk  L.  Yerkes  ^  Rebekah  Puterbauglf 
Air  Force  Research  Laboratory,  WPAFB  OH  45433 
Scott  K.  Thomas§ 

Wright  State  University,  Dayton,  OH  45435 
John  McQuillen^ 

NASA  Glenn  Research  Center,  Cleveland,  OH  44135 


This  paper  discusses  the  effects  of  a  variable-gravity  environment  on  the 
performance  of  a  subcooled  partially-confined  spray.  An  experiment,  consisting  of  a  test 
chamber,  the  associated  flow  loops,  and  instrumentation,  was  fabricated  to  be  flown  on  the 
NASA  Reduced-Gravity  Testing  Platform.  This  modified  KC-135  aircraft  was  flown 
following  a  parabolic  flight  path  to  provide  various  gravity  levels  of  approximately  0.16, 

0.35, 1.0,  and  1.8  g.  The  spray  chamber  contains  two  opposing  nozzles  spraying  on  Thick 
Film  Resistor  (TFR)  heaters,  which  are  mounted  on  insulating  glass  pedestals.  Only  the 
upward  facing  heater  was  used  during  this  testing.  Thermocouples  under  the  heater  in  the 
glass  pedestal  and  in  the  cooling  fluid  around  the  heater  are  used  to  determine  the  surface 
temperature  and  the  amount  of  heat  being  carried  away  by  the  cooling  fluid.  The  glass 
pedestals  are  surrounded  by  an  annular  sump  system,  which  is  used  to  collect  and  remove 
the  cooling  fluid  from  the  test  chamber.  The  fluid  used  for  this  testing  is  FC-72,  a 
refrigerant  that  is  non-toxic,  non-flammable,  and  non-reactive.  Due  to  its  dielectric  nature, 
FC-72  can  also  be  sprayed  directly  on  the  electric  thick  film  heaters.  For  this  set  of  testing, 
the  parametric  ranges  are  as  follows:  The  non-dimensional  heat  input  was  varied  over  the 
range  35  <  (\-f)GA  <  190,  which  corresponds  to  a  heat  flux  to  the  spray  of  20  <  (1  -f)q”  <  70 
W/cm2.  The  non-dimensional  grouping  (Fr1/2Ga)1/2  was  varied  over  the  range  35  < 
(Fr1/2Ga)12  <  70,  by  the  varying  acceleration  field  on  board  the  aircraft  and  by  varying  the 
flow  rates  of  the  cooling  fluid.  The  working  fluid  is  subcooled,  with  the  subcooling 
temperature  ranging  from  24°  C  <  Tsc  <  30  °C.  The  non-dimensional  grouping  (Fr1/2Ga)12  is 
a  ratio  of  inertia  and  acceleration  forces  to  viscous  forces.  For  this  data,  which  comes  from 
one  flight  week  consisting  of  three  flights,  the  temperature  difference  AT  =  Ts-Tsat  (surface 
temperature  -  saturation  temperature)  is  generally  found  to  decrease  with  decreasing  values 
of  (Fr1/2Ga)12  while  the  heat  input  remains  constant.  Similarly,  from  the  data  presented 
here,  it  appears  that  higher  values  for  subcooling  may  enhance  this  acceleration  effect.  Also, 
changing  the  flow  rate  may  affect  how  much  the  acceleration  level  affects  the  2-phase  heat 
transfer  coefficient  (/f2-<p),  with  the  higher  flow  rate  showing  less  of  a  change  in  /f2(p  as  the 
acceleration  level  changes. 


I.  Nomenclature 


a 

A 

a9b9c 

b 


=  acceleration  level,  g 
=  heater  area,  m2 

=  constants  to  be  determined  in  surface  temperature  formulation 
=  radius  of  the  glass  heater  post  assembly,  m 
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-  droplet  diameter,  m 

f 

=  heater  conduction  loss  fraction 

Fr 

=  Froude  number,  v2/aDv 

g 

=  volumetric  heat  generation,  W/m3,  Q/nb2Hh{v 

G 

=  gb2/((Tsat  -  TX'W an)  khtr) 

GA 

=  non-dimensional  heat  input,  q/( n  b  Arhtr  (Tsat  -  T0 0?wan)) 

Ga 

=  Galileo  number,  aDp p2/p2 

h 

=  convective  heat  transfer  coefficient,  W/(m2-K) 

^top 

=  convective  heat  transfer  coefficient  at  surface  of  heater,  W/(m2-K) 

^l-(p 

=  single-phase  convective  heat  transfer  coefficient,  W/(m2-K) 

^2-cp 

=  two-phase  convective  heat  transfer  coefficient,  W/(m2-K) 

HCOY 

=  thickness  of  heater  cover,  m 

Hhtr 

=  thickness  of  heater  layer,  m 

TTsub 

=  thickness  of  heater  substrate,  m 

k 

=  thermal  conductivity,  W/(m-K) 

kcov 

=  thermal  conductivity  of  heater  cover,  W/(m-K) 

kf\u\d 

=  thermal  conductivity  of  FC-72  coolant,  W/(m-K) 

khtr 

=  thermal  conductivity  of  heater  layer,  W/(m-K) 

^sub 

=  thermal  conductivity  of  heater  substrate,  W/(m-K) 

K 

—  ^fluid/^htr 

Nu 

=  Nusselt  number,  htop  b/kf\md 

T^sat 

=  saturation  pressure,  N/m2 

e 

=  heater  input  power,  W 

Qcov 

=  power  through  heater  cover,  W 

Qsub 

=  power  through  heater  substrate,  W 

q” 

=  heater  input  heat  flux,  W/cm2 

q  sen 

=  maximum  sensible  heat  flux,  W/cm2 

q  i-cp 

=  single-phase  heat  flux,  W/cm2 

q”  2-cp 

-  two-phase  heat  flux,  W/cm2 

T 

=  temperature,  °C 

Tint 

=  interface  temperature,  °C 

Th 

=  temperature  between  substrate  and  heater,  °C 

Tf 

=  free  stream  fluid  temperature,  °C 

Tnoz 

=  nozzle  inlet  temperature,  °C 

Ts 

=  heater  surface  temperature,  °C 

^sat 

=  saturation  temperature,  °C 

^sc 

=  subcooling  temperature,  Tsat-Tnoz ,  °C 

^sub 

=  substrate  temperature,  °C 

Tt 

=  temperature  between  heater  and  cover,  °C 

^GO,tOp 

=  local  freestream  temperature  of  liquid  flowing  over  heater  surface, 

Tco,wall 

=  local  freestream  temperature  of  liquid  flowing  along  post  wall,  °C 

V 

=  velocity,  m/s 

F 

=  volumetric  flow  rate,  m3/s 

We 

=  Weber  Number,  p v2Dv/  a 

X 

Greek  symbols 

=  position  variable  used  in  surface  temperature  formulation 

AT 

=  temperature  difference,  ( Ts  -  rsat),  °C 

6 

=  non-dimensional  temperature,  ( T  -  wan)  /  (rsat  -  T^waii) 

$int 

=  (Tint  -  T^wall)  /  (Tsat  “  T^wall) 

=  (Ts  -  TaO'W all)  /  (rsat  -  Foo  wall) 

^sat 

=  (Tsat  -  TTo  wall)  /  (Tsat  “  T^wall) 

^oo, top 

=  (Too, top  "  TV^wall)  /  (Tsat  -  T^wall) 

P 

=  density,  kg/m3 

a 

=  surface  tension,  kg/s2 
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II.  Introduction 

The  thermal  management  of  high-power,  high- flux  devices,  such  as  those  used  in  power  conditioning, 
power  electronics,  and  directed  energy  applications,  is  becoming  increasingly  difficult.  Even  with  advances  in 
materials  and  efficiencies,  these  devices  are  requiring  the  dissipation  of  ever  increasing  amounts  of  thermal  energy. 
Many  of  these  devices  also  have  low  thermal  masses  and  have  isothermality  requirements  of  ±  1  -  2  K.  Asa  result, 
robust  thermal  management  approaches  need  to  be  developed.  Various  methods  have  been  considered,  with  just  a 
few  of  these  being  two-phase  flow  and  heat  transfer1,2  and  pool  boiling.3”5  One  other  potential  solution  to  the 
thermal  management  of  these  devices  has  been  spray  cooling. 

While  spray  cooling  is  a  method  that  has  been  used  often  in  the  past,  there  has  been  little  research  done  on  the 
effects  of  a  variable  gravity  environment  on  the  effectiveness  of  spray  cooling.  Yoshida  et  al.6  and  Kato  et  al.7  both 
present  some  results  on  this  topic.  Some  of  the  more  recent  work,  however,  was  done  by  Baysinger  et  al.8  These 
researchers  presented  the  design  and  preliminary  testing  of  a  variable-gravity  spray  cooling  experiment.  The 
primary  purpose  of  this  study  was  to  investigate  heat  transfer  and  fluid  management  issues  for  a  continuous-flow, 
closed-loop  spray  cooling  system  subjected  to  a  variable-gravity  environment.  Tests  were  conducted  using  ITO 
(Indium  Tin  Oxide)  heaters  onboard  the  NASA  KC-135  Reduced-Gravity  Research  Aircraft,  which  provided  the 
variable-gravity  environment  by  following  a  parabolic  flight  trajectory.  This  preliminary  testing  provided  much 
information  on  the  fluid  management  aspects  of  variable  gravity  testing.  More  information  on  this  test  rig  can  be 
found  in  Baysinger.9 

Building  on  this  work  by  Baysinger  et  al.8  was  the  research  presented  by  Yerkes  et  al.10  This  work  was  again 
performed  on  the  NASA  KC-135  with  ITO  heaters,  but  with  higher  heat  loads  than  the  experiments  that  were 
presented  previously.  The  non-dimensional  heat  input  for  these  tests  was  varied  from  (\-f)GA  =  30  to  110,  and  the 
non-dimensional  grouping  (Fr1/2Ga)12  was  varied  from  20  to  66.  During  this  testing,  the  non-dimensional 
temperature  difference  was  found  to  decrease  with  decreasing  (Fr1/2Ga)12  for  a  given  heat  load,  suggesting  that  the 
Nusselt  number  was  actually  increased  by  decreasing  (Fr1/2Ga)1/2.  However,  this  testing  was  primarily  in  the  single¬ 
phase  region.  It  was  also  determined  in  this  testing  that  the  ITO  heaters  were  prone  to  failure  at  the  higher  heat 
loads,  due  to  high  current  densities  in  the  thin  film  of  resistive  ITO  where  the  heating  actually  takes  place. 

The  objective  of  the  present  investigation  is  to  continue  the  work  presented  by  Baysinger  et  al.8  and  Yerkes  et 
al.10,  and  to  investigate  the  effects  of  variable-gravity  on  the  cooling  performance  of  a  partially-confined  FC-72 
spray.  The  experiment  was  again  carried  out  using  the  NASA  KC-135  to  supply  the  variable  gravity  environment, 
but  the  ITO  heaters  were  replaced  by  thick  film  resistor  (TFR)  heaters,  which  are  capable  of  handling  the  higher  heat 
loads.  The  cooling  performance  was  tested  by  varying  the  non-dimensional  heat  input  35  <  (l-f)GA  <190  for  the 
non-dimensional  grouping  35  <  (Fr  Ga)  <  70.  This  corresponds  to  a  heat  flux  range  of  20  <  <  70  W/cm  . 

The  range  of  subcooling  for  these  tests  was  24°C  <  Tsc  <  30°C. 

III.  Experimental  Design 

The  experimental  apparatus  used  in  the  current  research  is  very  similar  to  that  described  by  Baysinger  et  al.8. 
This  experiment  consists  of  a  test  rig,  which  houses  the  fluid  management  and  data  acquisition  systems,  and  a  test 
chamber,  which  houses  two  opposing  spray  nozzles,  the  heater  pedestal  assemblies,  and  the  fluid  control  /  removal 
sump  devices  around  the  pedestals.  Figure  la)  shows  the  entire  rig,  while  figure  lb)  shows  the  spray  chamber. 
Only  the  upward  facing  heater  was  used  during  this  testing. 


Figure  1.  Experimental  test  rig  showing  the:  ( a )  complete  package  and  (b)  test  chamber. 
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The  test  rig  contained  three  flow  loops:  a  water  loop,  a  spray  nozzle  loop,  and  a  chamber  drain  loop.  The  water 
loop  on  this  rig  is  used  only  for  control  of  the  temperature  and  pressure  environment  inside  the  chamber.  The  drain 
loop  was  used  for  draining  the  chamber,  if  the  sumps  were  unable  to  remove  sufficient  fluid.  The  nozzle  loop  was 
used  to  actually  deliver  the 
working  fluid,  FC-72, 
through  the  spray  nozzles 
and  onto  the  heater  surface. 

FC-72  was  chosen  as  the 
working  fluid  for  various 
reasons.  It  is  non-toxic, 
non-flammable,  non- 
corrosive,  and,  as  a 
dielectric,  it  can  be  sprayed 
directly  on  the  electrical 

component  to  be  cooled  (in  Figure  2:  Pedestal  photos  showing:  a)  close-up  of  heater  region  and  b)  entire 

this  case,  the  TFR  heater).  pedestal. 

The  heaters  that  are  used  in  this  set  of  testing  are  similar  to  the  heater  /  pedestal  assembly  described  in  Yerkes  et 
al.10.  These  cylindrical  pedestals  consist  of  several  layers.  The  first  three  layers  comprise  the  TFR  heater.  The 
layers  of  the  TFR  heater  were  investigated  by  Glaspel11.  Table  1  summarizes  the  sizes  and  thermal  conductivities  of 
these  heater  layers.  These  heater  layers  sit  on  top  of  two  thin  wafers  of  glass  which  are  each  about  0.001  m  thick. 
Finally,  these  glass  layers  sit  on  top  of  a  glass  pedestal,  which  supports  and  insulates  the  heater,  in  an  attempt  to 
minimize  heater  loss.  A  small  (0.0254  cm  diameter)  Type  E  thermocouple  is  embedded  between  the  heater  substrate 
and  the  first  glass  layer  of  the  pedestal.  Several  others  of  these  thermocouples  are  embedded  between  the  glass 
layers.  Together,  these  thermocouples  allow  for  calculations  of  how  much  heat  is  being  lost  through  conduction 
down  the  glass  pedestal,  and  what  the  surface  temperature  is  at  the  top  of  the  glass  cover  over  the  heater  film.  From 
the  results  of  previous  experiments,  the  percentage  of  heat  being  lost  down  the  pedestal  is  estimated  to  be /=  1.5% 
of  the  total  heat  supplied  by  the  heater.  The  pedestal  was  surrounded  by  an  annular  sump  structure  fitted  with  a 

containment  cap. 

This  sump  and  cap  served  both  as  fluid  containment  and 
as  a  means  of  removing  the  fluid  from  the  chamber.  Two 
thermocouples  were  positioned  inside  this  sump  annulus  to  try 
to  monitor  the  temperature  of  the  fluid  that  was  coming  off  the 
heater.  One  thermocouple,  the  upper  film  temperature,  was 
placed  just  off  the  surface  of  the  heater,  and  a  second,  the 
lower  film  temperature,  was  about  1.2  cm  down  from  the 
upper  film  temperature,  near  the  side  of  the  glass  pedestal. 
These  temperatures,  along  with  the  nozzle  inlet  temperature, 
allow  the  calculations  of  the  average  film  temperature  that  can 
be  used  with  the  heater  surface  temperature  to  determine  the  heat  transfer  coefficient  of  the  spray  cooling.  The 
nozzle  inlet  temperature,  and  therefore  the  subcooling  of  the  spray,  can  be  controlled  by  a  preheater  on  the  flow  loop 
leading  up  to  the  nozzle  inlet.  This  preheater  consisted  of  a  length  of  copper  tubing  wrapped  with  a  heater  tape, 
controlled  by  a  PID  temperature  controller.  This  preheater  could  be  used  to  bring  the  temperature  of  the  working 
fluid  coming  out  of  the  nozzle  closer  to  the  saturation  temperature  in  the  chamber.  This  saturation  temperature  was 
dictated  by  the  pressure  in  the  chamber,  which 
could  be  controlled  by  the  temperature  of  the 
water  flowing  around  the  chamber  in  the  water 
loop. 

The  thermocouples  in  and  around  the 
pedestal  in  the  test  chamber  are  all  referenced 
through  a  dry-well  ice  point  to  improve  the 
stability  and  accuracy  of  the  temperature 
readings.  Additionally,  the  pedestal 

thermocouples  are  calibrated  using  an  oil  bath  / 

RTD  calibration  setup.  The  acceleration  of  the 
variable  gravity  environment  is  measured  by  a  tri-axis  accelerometer  mounted  directly  to  the  test  chamber.  All  of 
the  data  that  are  collected  from  the  thermocouples,  the  flow  meters,  and  the  pressure  transducers  on  the  flow  loops 
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Table  1:  Heater  Layer  Dimensions  and  Conductivities 


Layer 

Thickness  (m) 

Thermal 
Conductivity 
(W  /  m-K) 

Ceramic  Substrate 

Hsuh  =  0.000634 

ksub  =  27.0 

Heater 

Hhtr  =  0.000008 

h tr=  1-04 

Glass  Cover 

Hcoy  =  0.000040 

kCOv  =  1  -04 

Figure  3:  Pedestal  in  sump  and  cap 


are  routed  through  a  data  acquisition  system  to  a  laptop.  The  experiment  is  controlled  using  a  control  panel  at  the 
top  of  the  test  rig,  next  to  the  spray  chamber. 

The  heater  power  was  measured  using  two  different  voltage  measurements.  The  first  was  across  the  heater  film 
itself,  as  close  to  the  heater  as  possible  to  minimize  voltage  error  from  the  heaters  leads.  The  second  voltage 
measurement  was  across  a  precision  resistor.  This  measurement  allowed  for  the  calculation  of  current  (current  = 
voltage  /  resistance).  This  current  could  then  be  multiplied  by  the  voltage  across  the  heater  to  give  the  heater  power. 

The  data  presented  here  was  collected  during  one  flight  week  on  the  NASA  KC-135,  which  consisted  of  three 
flights  of  about  forty  parabolas  each.  Each  parabola  provided  approximately  35  seconds  of  elevated  gravity, 
followed  by  about  25  seconds  of  reduced  gravity.  To  collect  the  data,  the  flow  rate  and  level  of  subcooling  were  set 
at  the  beginning  of  the  flight,  and  the  heater  power  input  was  varied  to  the  desired  non-dimensional  heat  inputs.  The 
heat  level  was  maintained  for  at  least  one  full  parabola  (elevated  and  reduced  gravity),  and  then  was  adjusted  to  the 
next  desired  level.  The  set  flow  rate  and  the  varying  accelerations  gave  the  varying  levels  for  the  non-dimensional 
grouping  (Fr1/2Ga)1/2.  Table  2  shows  the  ten  various  cases  of  flow  rate,  subcooling,  and  the  non-dimensional 
grouping  (Fr1/2Ga)1/2  for  these  three  flights. 


Table  2:  Cases  for  the  heat  flux  versus  temperature  difference  plots 


Case 

TS.,CC ) 

TK  (°C) 

Acceleration 

(g> 

(Fr^'Ga)1^ 

Volumetric 
Flow  Rate 

V  *  10*  u+%) 

Weber 

Number 

Flight 

l 

Case  1 

53.7  ±1.0 

25.2  +  1.1 

0.15  ±0.02 

38.26+  1.26 

8.70  +  0  11 

796  +  17 

Case  2 

53.7  ±0.9 

34.9  ±  1.0 

1.04  +  0  15 

61 .57  ±  2.53 

8  63  +  0  14 

787  ±19 

Case  3 

53.5  ±1.1 

24  6+  12 

1.78  +  0  04 

71  26  +  1-00 

8.73  +  0  1  ] 

$04  +  18 

Flight 

Case  4 

60.0  ±0.3 

31.3  ±0.3 

0.37  +  0  04 

47,25  ±  118 

8  62  +  0  09 

782  +  14 

Case  5 

59.$  ±0.4 

30.3  ±0.5 

1  06  +  0  15 

61 .84  ±  2.34 

8  71  ±005 

795  +  7 

Case  6 

59.8  ±0.5 

30.8  +  0.' 

1  "8  +  0  06 

70.64  ±09 

8.60  +  0.07 

"81  +12 

- - - 

Flight 

3 

Case  7 

59.4  +  0.4 

30.4  +  0." 

0.16  +  0  03 

34  89+  1  48 

6  54  +  00" 

522  +  7 

Case  8 

59.7  ±0.2 

30.8+0.4 

0.36  +  0.04 

42.87  +  1  12 

6.57  +  0.10 

525  +  7 

Case  9 

58.8  ±0.3 

29.5  ±0.5 

1.06  +  0.15 

55.91  +  2.03 

6.60  +  0.08 

527  +  8 

Case  10 

59.3  ±  0  3 

29.7  ±05 

1.76  +  0  05 

63.54  ±  0.57 

6  44  +  0.08 

514  +  “ 

IV.  Mathematical  Formulation 

A.  Surface  temperature  calculations: 

The  ITO  (Indium  Tin  Oxide)  heaters  that  were  originally  used  on  the  pedestals  were  not  capable  of 
reaching  CHF  (critical  heat  flux)  for  the  flow  rates  that  were  being  tested.  If  the  power  was  increased  to  try  to  reach 
CHF,  the  heaters  failed  due  to  the  high  current 
density  in  the  thin  film  of  the  ITO  heater.  These 
heaters  were  replaced  by  thick  film  resistor 
(TFR)  heaters.  This  heater  consisted  of  a 
ceramic  substrate,  a  thick  film  resistive  element, 
and  a  glass  cover  plate.  The  thicknesses  and 
thermal  conductivities  of  these  materials  are  in 
Table  1. 

Figure  4  shows  the  geometry  of  the  top 
of  the  heater  pedestal  (not  to  scale).  The  heat 
transfer  through  the  substrate  and  cover  are  also 
given  on  this  schematic,  where  Q  is  the  power 
produced  by  the  resistive  heater  and  /  is  the 
fraction  of  heat  lost  down  the  pedestal  (assumed 
to  be  1.5%,  as  in  previous  papers).  The  dashed 
line  on  this  plot  represents  the  assumed 
temperature  distribution  through  this  TFR  heater. 


Figure  4:  Heater  geometry  with  heat  flow  and 
temperature  distribution 


This  is  based  on  one-dimensional  heat  transfer  through  the  heater.  For  the  purposes  of  reducing  the  data,  the  heater 
power,  Q ,  and  the  temperature  at  the  interface  between  the  substrate  and  the  glass  pedestal,  rint,  are  known.  The 
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temperature  profiles  in  the  substrate  and  cover  regions  are  assumed  to  follow  one-dimensional  conduction,  and  the 
temperatures  at  the  interfaces  are  the  only  relevant  temperatures  in  this  case.  The  following  equations  show  the 
equations  for  these  temperatures. 

fT~T,t  \ 


<2sub  =-K^A 


dT  _Tb-Tmt 


a 


sub 


0COV  =  ~kcovA 


#sub 

J  dx 

#sub 

Ub 

Ts~Tt) 

dT 

— s — 

1 

Uh 

1 

Qcov 

ZJ 

11  COV  J 

— y  -  — 

dx 

^cov 

Ko,A 

(1) 


7b,  the  temperature  between  the  substrate  and  the  heater,  can  be  calculated  as: 

fQH 


Tb=- 


Qsub  ^^sub  _|_  j1  _  _ J  '^1  ±  sub  _|_ 


Akt 


mt  Ak 

sub  'Nub 


(2) 


The  temperature  distribution  within  the  heater  can  be  assumed  to  be  a  parabolic  function  in  terms  of  the  X  shown 
on  the  schematic.  This  parabolic  temperature  has  to  match  the  temperatures  and  slopes  set  by  the  heat  transfer 
through  the  substrate  and  cover: 

T (x)  =  Ax 2  +  Bx  +  C ;  0  <  x  <  Hh 


■  htr 


x  =  0,  T(0)  =  Tb 

-=o,  a* 

dx  Ak. 


(3) 


X=Hte> 


dT 


sub 

a 


dx  Akcm. 
x=HhtI,  T(Hht[)  =  Tt 

Using  the  first  three  boundary  conditions  (BC)  to  solve  for  the  unknown  constants  v4,5,  C  and  then  substituting 
these  back  into  the  distribution  gives: 


T(x)  = 


1 


2  AH, 


htr 


^sub  a 


^sub 
V  ^sub 


X 


®mh-X  +  T 


COY  J 


Ak. 


(4) 


sub 


Now,  substitute  in  X  =  Hh{r  (fourth  BC)  to  get  the  value  for  Tt,  the  temperature  at  the  interface  between  the  heater 
layer  and  the  glass  cover: 


T  =- 


1 


2  AH 


htr 


Qsuh  Qo 


^sub 

V  Ub 


(Hj2- 


Gsub^h 


'sub  htr  _|_  y 


cov  J 


Ak. 


sub 


The  temperature  at  the  surface  of  the  glass  cover  can  be  given  as: 

=T,  (1 
1 1 


T  =T. 


(5) 


(6) 


(7) 


Ak»,  AKov 

Now,  substituting  in  the  value  of  Tt  given  above  and  simplifying  gives  the  surface  temperature  as: 

a  -rmw 

Akcov 

This  is  the  surface  temperature  that  is  used  for  the  rest  of  the  data  reduction  calculations. 

B.  Two-Phase  Heat  Transfer  Coefficient  Definition 

The  standard  definition  of  a  single-phase  (1-cp)  heat  transfer  coefficient  is  derived  from  the  equation  for 
convective  heat  transfer,  which  is: 


T  =Tb  +  ®HhtI  ■ 

f 

1 

1 

+ 

(1-/111 

( 

s  b  A 

k  sub 

2 

LUb 

*cov  JJ 

~(p  1  ~(p  ' 

This  leads  to  a  heat  transfer  coefficient  that  can  be  defined  as 


(8) 
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This  equation  is  used  to  calculate  the  h\.^  for  a  single  point  of  the  data.  Another  way  of  looking  at  this  is  to 
plot  q ”  versus  ( Ts-Tf ).  The  average  heat  transfer  coefficient  can  then  be  viewed  as  being  the  slope  of  the  q ”  versus 
(Ts-Tf)  curve.  This  is  all  valid  for  a  1-  cp  system.  However,  for  this  testing,  and  due  to  the  subcooling,  the  heater 
surface  temperature  was  initially  below  the  saturation  point.  With  increasing  heat  input,  the  surface  temperature 
increased  until  it  was  above  the  saturation  point.  For  a  surface  temperature  above  the  saturation  point,  it  was 
assumed  that  there  is  vapor  being  generated,  resulting  in  a  two-phase  (2-cp)  system. 

The  convention  for  a  2-cp  system  is  to  plot  q ”  versus  ( Ts-Tsat ).  For  this  convention,  the  two  phase  heat 
transfer  coefficient,  h2.(?,  can  then  be  viewed  as  the  slope  of  q ”  versus  AT.  Also,  this  heat  transfer  coefficient  can  be 
determined  directly  using  equation  (12)  which  accounts  for  the  sensible  heating  occurring  at  Ts  =  Tsat.  Being  able  to 
calculate  the  h2_ 9  value  directly  from  the  data,  without  first  plotting  it  and  finding  a  curve  fit  and  slope,  can  be 
beneficial  under  certain  situations. 


Figure  5:  Heat  flux  versus  temperature  difference  plots  for  a)  single-phase  and  b)  two-phase 


Figure  5  b)  shows  the  heat  transfer  curve  with  heat  flux  plotted  against  the  temperature  difference  AT  =  ( Ts - 
rsat)  .  For  this  formulation,  it  was  assumed  that  the  region  of  AT  <  0  is  single-phase,  and  that  the  region  of  AT  >  0  is 
two-phase.  Thus,  the  1-  cp  portion  of  this  curve  is  basically  the  q ”  versus  (Ts-Tf)  plot,  figure  5  a),  but  transformed  to 
a  q ”  versus  A T=  Ts  -  Tsat  coordinate  system.  The  heat  flux  at  the  point  where  the  curve  crosses  the  q ”  axis  (A T=  0) 
is  assumed  to  be  the  maximum  sensible  heat  flux  that  the  1-  cp  fluid  can  achieve  (here,  called  q'\Q n),  prior  to  the 
onset  of  vapor  generation.  The  total  heat  flux  can  be  calculated  as: 

^'=V,(7;-7»|  +  h2_<p(.Ts-Tsat)^q''=q''sen+h2_<p(Ts-Tsat)  (10) 

us  ~1sat 

where  it  has  been  assumed  that 


Kr{T,~Tr) 


q 


(ii) 


This  assumption  is  valid  if  the  fluid  temperature  and  the  saturation  temperature  remain  constant.  This  is  the 
case  in  the  current  research.  Now,  equation  (10)  can  then  be  solved  for  h2< p  to  give: 


h  q  q  sen 

2-cp 

s  ±  sat  ) 


(12) 


This  equation  can  then  be  used  to  calculate  the  two-phase  heat  transfer  coefficient  on  a  point  by  point  basis, 
much  like  equation  (9)  calculated  the  single-phase  heat  transfer  coefficient  at  a  given  point.  The  slope  of  the  two- 
phase  portion  of  the  q”  versus  AT  curve  can  give  the  average  value  for  the  h2_ 9  over  the  entire  two-phase  test  case. 


IV.  Results  and  Discussion 

Experimental  runs  were  conducted  over  the  course  of  three  flight  tests  as  shown  in  Table  2.  The  FC-72 
spray  was  saturated  with  air,  approximately  30-40%  by  volume,  based  on  a  partial  pressure  calculation,  to  minimize 
maintenance  of  the  thermal  system  in  the  event  of  a  needed  repair.  However,  operating  the  spray  system  with  air 
saturated  FC-72  resulted  in  an  elevated  saturation  temperature  that  limited  the  maximum  temperature  to  below  onset 
of  CHF.  Table  2  describes  the  ten  different  parameter  configurations  that  were  tested  during  this  set  of  three  flights. 
For  each  flight,  a  given  Tsah  Tsc ,  and  V  were  set  and  maintained  at  constant  levels  throughout  the  duration  of  the 
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flight.  The  non-dimensional  quantity  (Fr1/2Ga)1/2  was  varied  due  to  the  varying  acceleration  level  on  the  aircraft. 
Figures  6-8  show  the  heat  flux  input  to  the  spray  versus  the  temperature  difference,  AT  =  Ts  -  rsat,  for  these  ten  cases. 
These  plots  are  presented  in  both  dimensional  and  non-dimensional  form. 

A.  Heat  Flux  Versus  Temperature  Difference 

Figure  6  a)  shows  the  dimensional  (1  -f)q”  versus  A T=  Ts-Tsat  plot  for  flight  1,  while  figure  6  b)  gives  the 
non-dimensional  (l-f)GA  versus  0s-0sat  plot  for  this  same  flight.  In  figure  6  a),  the  AT  increases  with  increasing  heat 
flux,  as  would  be  expected.  The  left  hand  side  of  the  plot  (with  the  (1  -f)q”  axis,  AT  =  0,  as  the  separating  line)  is 
taken  to  be  the  single-phase  region,  and  the  right  hand  side  of  the  plot  is  taken  to  be  two-phase.  The  data  is  plotted 
and  linear  trend  lines  are  applied  to  the  l-(p  and  2-cp  regions  separately.  There  is  a  pronounced  shift  in  the  data  for 
the  three  different  cases.  For  case  1,  with  the  lowest  value  of  (Fr1/2Ga)1/2,  the  AT  value  for  any  given  heat  flux  is 
significantly  less  than  the  AT  values  for  both  case  2  and  case  3  at  the  same  heat  flux.  The  maximum  difference 
between  the  values  of  both  V  and  We  for  these  three  cases  is  less  than  3%,  and  the  Tsc  and  Tsat  values  also  change 
very  little,  on  the  order  of  2%.  The  acceleration  and  the  (Fr1/2Ga)12  value,  on  the  other  hand,  change  drastically 
between  the  three  cases,  and  this  is  the  only  real  difference  between  the  parameters  for  cases  1,  2,  and  3.  There  is 
some  shift  between  cases  2  and  3,  but  it  is  not  as  pronounced  as  between  cases  1  and  2.  This  means  that  between  a 
=  0.15  g  and  a  =  1.04  g,  there  is  a  large  change,  but  for  between  a  =  1.04  g  and  a  =  1.78  g  there  is  much  less 
change.  Another  point  to  note  is  that,  for  this  figure,  there  is  very  little  difference  between  the  slopes  of  the  1-cp  and 
2-cp  regions.  The  data  appear  to  be  very  linear  throughout  the  entire  range  of  AT.  Also,  the  curves  are  relatively 
parallel.  The  curves  diverge  a  small  amount  as  the  AT  increases,  but  the  slopes  of  the  three  cases  are  approximately 
the  same.  Many  of  the  same  trends  can  be  seen  in  figure  6  b).  The  non-dimensional  temperature  difference  0s-0sat 
tends  to  increase  as  the  non-dimensional  heat  flux  ( \-f)GA  increases.  There  is  still  a  difference  between  cases  1  and 
2,  showing  that  the  lower  acceleration  of  case  1  drops  the  value  of  0s-0sat  at  a  given  ( \-f)GA  when  compared  to  case 
2.  However,  there  is  a  slightly  more  pronounced  difference  between  case  2  and  case  3  when  the  non-dimensional 
data  is  observed,  as  compared  to  the  dimensional  data  of  figure  6  a).  The  slopes  in  figure  6  b)  behave  very  similarly 
to  those  observed  in  6  a),  with  the  curves  being  relatively  linear  and  parallel. 

The  plots  for  flight  2  are  given  in  figure  7,  with  figure  7  a)  showing  the  dimensional  data  and  figure  7  b) 
showing  the  non-dimensional  data.  Again,  the  plots  show  both  the  data  and  linear  trend  lines,  split  up  between  the 
1-cp  and  2-cp  regions.  In  this  case,  there  is  some  difference  between  the  slopes  of  the  1-cp  and  2-cp  regions,  though  it 
is  still  not  a  drastic  change.  The  curves  still  appear  to  be  generally  linear  throughout  the  entire  range  of  AT.  Also, 
the  curves  are  still  relatively  parallel.  When  viewing  the  three  different  cases  of  figure  7,  there  is  again  a  shift  in  the 
AT  and  #s-0sat  data.  For  case  4,  which  is  at  a  lower  acceleration  (a  =  0.37  g),  the  temperature  difference  is  smaller  at 
a  given  heat  flux  when  compared  to  case  5.  However,  in  both  the  dimensional  and  non-dimensional  plots,  there  is 
more  of  a  difference  between  case  5  (a  =  1.06  g)  and  case  6  (a  =  1.78  g)  than  between  case  2  and  case  3  in  figure  6. 

Figure  8  shows  the  data  from  flight  3.  There  is  little  difference  between  the  two  reduced-acceleration  cases 
(case  7,  at  a  =  0.16  g,  and  case  8,  at  a  =  0.36  g).  The  same  can  be  said  for  case  9  (a  =  1.06  g)  and  case  10  (a  =  1.76 
g).  There  is  still  the  familiar  shift  between  cases  7  and  8  and  cases  9  and  10.  In  figure  8,  however,  the  curves  do  not 
appear  to  be  as  parallel  as  in  figures  6  and  7.  There  is  a  good  deal  of  divergence  between  cases  7  and  8  and  cases  9 
and  10  as  AT  increases.  This  suggests  that  there  is  a  noticeable  increase  in  the  convective  heat  transfer  coefficient  as 
the  acceleration  level  decreases. 

In  all  of  these  figures,  the  same  general  trend  is  observed.  For  lower  values  of  the  non-dimensional 
quantity  (Fr1/2Ga)1/2  (i.e.,  lower  accelerations),  the  temperature  difference  tends  to  be  lower  for  a  given  heat  input. 
One  initial  thought  was  that  it  may  have  been  a  flow  rate  issue,  where  the  flow  rate  was  increasing  at  the  lower 
acceleration  levels  (lower  (Fr1/2Ga)1/2).  However,  looking  at  table  2,  the  flow  rate  varies  by  less  than  3%  for  the 
different  cases  across  each  flight.  For  instance,  the  flow  rates  for  the  three  different  cases  on  flight  one  were  all 
within  about  0.1  m3/s  of  each  other.  This  suggests  that  the  lowered  surface  temperature  trend  must  be  due  to 
something  besides  a  change  in  flow  rate.  The  saturation  temperature  and  the  nozzle  inlet  temperature  were  both  kept 
fairly  constant  over  the  duration  of  each  flight  (resulting  in  a  constant  level  of  subcooling).  It  is  concluded  that  the 
acceleration  is  a  major  driving  factor  behind  the  change  in  behavior  of  the  surface  temperature  (since  Tsat  remains 
constant). 

Looking  at  figures  6  and  7  (flights  1  and  2)  can  give  insight  into  the  effect  that  varying  the  subcooling  can 
have  on  the  acceleration  effects  for  this  system.  The  volumetric  flow  rate  for  both  of  these  flights  was 
approximately  8  m3/s,  but  the  subcooling  for  cases  4,  5,  and  6  was  higher  than  that  for  1,  2,  and  3.  On  figure  7,  the 
curves  are  spread  further  apart  than  on  figure  6,  as  mentioned  above.  In  other  words,  the  reduction  in  the 
temperature  difference  is  more  pronounced  in  the  cases  with  the  higher  subcooling.  This  suggests  that  the  higher 
subcooling  actually  seems  to  enhance  the  observed  acceleration  effect.  Also,  there  is  more  of  a  slope  change 
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between  the  single-phase  and  two-phase  portions  of  the  curves  in  figure  7,  although  the  curves  are  still  relatively 
linear  throughout  the  entire  range  of  heat  fluxes  tested. 

Comparing  figures  7  and  8,  the  effect  that  varying  the  flow  rate  has  on  the  acceleration  effects  can  be 
examined.  In  both,  there  is  a  good  deal  of  separation  between  the  low-acceleration  curves  and  the  higher- 
acceleration  curves,  as  would  be  expected  with  the  higher  subcooling.  However,  in  figure  7,  the  curves  are  all 
relatively  parallel,  especially  in  the  2-phase  region.  There  is  some  slope  change  at  the  ATM)  point,  but  the  slopes  are 
basically  the  same  between  the  three  curves.  On  figure  8,  it  appears  that  the  lines  are  much  less  parallel,  especially 
on  the  two-phase  half  of  the  curve.  This  means  that,  in  flight  3,  there  is  a  higher  h2. 9  in  the  lower  accelerations, 
since  the  slope  of  these  curves  is  an  indication  of  the  convective  heat  transfer  coefficient.  One  explanation  for  this  is 
that  the  higher  flow  rate  may  tend  to  dampen  out  some  of  the  effect  that  the  acceleration  has  on  the  heat  transfer 
coefficient  in  the  two-phase  region.  In  all  three  figures,  the  heat  transfer  coefficient  (the  slope)  in  the  single-phase 
region  is  higher  for  the  lower  accelerations,  which  would  explain  why  there  is  such  a  temperature  difference 
between  the  lower-acceleration  curves  and  the  higher-acceleration  curves.  It  is  only  in  figure  8,  with  the  lower  flow 
rate,  that  this  higher  heat  transfer  coefficient  carries  over  into  the  two-phase  region. 


9 

American  Institute  of  Aeronautics  and  Astronautics 


(MGA  (!-/)</”  (W/cm2) 


A  r(°c) 


<>  case  u-rh 

“Case  2,1-Ph 
£k  Case  3,i-Ph 

*  Case  1,2 -Ph 
■  Case  2,2-Ph 

*  Case  3,2  Ph 
Linear  Kits 

"C  ase  LI -Hi 

Case  2J-Ph 
— Case  3,1-Ph 
Case  1,2-Ph 
“  "Case2,2-Ph 
Case  3,2  Ph 


-0.6  -0.4  -0.2  0.0  0.2  0.4  0.6  0.8  1.0  1.2 


0-0t 


sat 


Figure  6:  Heat  flux  to  the  spray  versus  temperature  difference,  Jsc  =  24.8°C,  Jsat  =  53.6°C, 
V  =  8.69*  10"6  m3/s,  showing  a)  dimensional  and  b)  non-dimensional  formulations 
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Figure  7:  Heat  flux  to  the  spray  versus  temperature  difference,  Tsc  =  30.7°C,  Tsat  =  59.8 
V=  8.64*  10'6  m3/s,  showing  a)  dimensional  and  b)  non-dimensional  formulations 
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Figure  8:  Heat  flux  to  the  spray  versus  temperature  difference,  Jsc  =  29.9°C,  Jsat  =  59. 
V  =  6.47*  10"6  m3/s,  showing  a)  dimensional  and  b)  non-dimensional  formulations 
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B.  Heat  Transfer  Coefficient  Comparison 

Table  3  shows  the  slopes  and  intercepts  for  the  various  linear  trend  lines  from  figures  6  a),  7  a),  and  8  a).  It 
also  shows  the  average  and  standard  deviation  of  the  calculated  h2< p  value,  and  the  percentage  difference  between  the 
2-cp  slope  and  the  average  calculated  h2< p  value.  For  each  case,  the  calculated  h2< p  value  was  calculated  using 
equation  (12),  with  the  g”se n  value  taken  as  the  average  of  the  1-cp  and  2-cp  intercept  value  for  that  case. 


Table  3:  Comparison  of  heat  transfer  coefficient  from  calculation  and  slope 


Case 

10000  *  1-cp 
Slope  From 
q "  vs.  AT 
(W/(m2-K)) 

1-cp 

Intercept 
From  q "  vs. 
AT  (W/cm2) 

10000  *  2-cp 
Slope  From 
q"v s.  AT 
(W/(m2-K)) 

2-cp 

Intercept 
From  q "  vs. 
AT  (W/cm2) 

Average 

Calculate 

d  h  7.<p 

(W/m2-K) 
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Deviation  of 
Calculated 

h  ?_(p 

Difference: 
Slope  and 
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h2-r  (%) 
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1 
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16300 

40.3 
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39.9 

14800 

2100 

7.5 
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14500 

36.5 

14700 

36.7 

14400 

1500 

o 

ci 

Case  3 

14500 

35.9 

15100 

35.3 

14900 

1600 

1.3 

Flight 

o 

Case  4 

17400 

51.3 

19100 

50.6 

18800 

3000 

1.6 

Case  5 

14100 

43.5 

20400 

43.9 

23200 

9300 

13.7 

Case  6 

13800 

40.7 

20100 

39.1 

18600 

2000 

7.5 

Flight 

3 

Case  7 

13900 

39.9 

23600 

37.7 

26200 

11200 

11.0 

Case  8 

16400 

42.1 

18900 

41.7 

18600 

1500 

1.6 

Case  9 

12800 

35.6 

17000 

34.9 

16400 

1900 

3.5 

Case  10 

13000 

36.1 
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35.4 

15000 

1400 

3.2 

From  table  3,  it  can  be  seen  that  the  2-phase  heat  transfer  coefficient  calculated  from  the  slope  of  the  q” 
versus  AT  plot  is  very  similar  to  the  average  h2_ 9  calculated  using  equation  (12).  The  maximum  percentage 
difference  between  these  two  values  on  any  given  case  is  13.8  %.  This  might  appear  to  be  a  fairly  large  difference, 
at  first  glance.  However,  there  is  a  fairly  large  standard  deviation  in  the  calculated  h2< p  values.  In  all  ten  of  the 
cases,  the  slope  number  falls  well  within  the  range  defined  by  the  average  calculated  h2_ 9  plus  or  minus  the  standard 
deviation  on  that  calculated  value.  The  large  standard  deviation  observed  in  the  average  calculated  h2_ 9  is  due  to  the 
general  variability  in  the  raw  data.  This  variability  in  the  raw  data  is  not  as  evident  in  the  slope,  as  it  is  an  average 
value  over  the  entire  two-phase  region,  whereas  the  calculated  h2_ 9  is  an  average  value  over  just  the  raw  points 
making  up  each  individual  data  point.  Since  the  slope  is  the  heat  transfer  coefficient,  this  indicates  that  the  method 
outlined  above  for  calculating  h2< p  from  the  data  is  valid. 

V.  Conclusions 

The  effects  of  variable  gravity  on  the  cooling  performance  of  a  partially-confined  spray  have  been  investigated. 
Examining  the  data  from  this  flight  week  shows  that,  as  the  acceleration  (and  thus  the  value  of  (Fr1/2Ga)1/2)  is 
decreased,  there  is  a  significant  downward  shift  in  the  AT  values,  for  a  given  (1  -f)q”  value.  This  same  trend  is 
observed  in  the  non-dimensional  temperature  difference  and  non-dimensional  heat  input.  At  this  time,  the  physics 
driving  the  observed  phenomena  are  not  clearly  understood.  Similarly,  from  the  data  presented  here,  it  appears  that 
higher  values  for  subcooling  may  enhance  this  acceleration  effect.  In  addition,  changing  the  flow  rate  may  play  a 
role  in  the  extent  to  which  the  acceleration  level  affects  the  2-phase  heat  transfer  coefficient.  The  higher  flow  rate 
shows  less  of  a  change  in  h2.^  as  the  acceleration  level  changes.  Finally,  experimental  data  were  obtained  using  a 
saturated  FC-72  liquid.  It  is  not  known  how  dissolved  air  in  the  FC-72  affects  the  heat  transfer  with  varying  gravity. 

In  addition,  the  effect  of  variable  gravity  on  the  critical  heat  flux  of  the  spray  cooling  was  not  studied  in  this  set 
of  tests.  The  interface  temperature  was  limited  by  a  temperature  cut-off  switch  to  try  to  protect  the  heater  from  over¬ 
temperature  situations,  and  to  try  to  keep  the  FC-72  from  dissociating  into  toxic  compounds.  For  future  testing,  the 
limit  on  the  interface  temperature  can  be  increased,  and  the  saturation  temperature  decreased  by  varying  the  chamber 
pressure,  to  allow  the  critical  heat  flux  to  be  studied  under  variable-gravity  conditions. 
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